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Abstract 

Fits to high energy data alone cannot cleanly discriminate between asymptotic In s and In 2 s behavior 
of total hadronic cross sections. We demonstrate that this is no longer true when we require that these 
amplitudes also describe, on average, low energy data dominated by resonances. We simultaneously fit 
real analytic amplitudes to high energy measurements of: f ) the ir + p and ir~p total cross sections and 
p- values (ratio of the real to the imaginary portion of the forward scattering amplitude), for yfs > 6 
GeV, while requiring that the asymptotic fits smoothly join the cr T + p and <r^_ p total cross sections at 
\fs~ =2.6 GeV — both in magnitude and slope , and 2) separately simultaneously fit the pp and pp total 
cross sections and p-values for *Js > 6 GeV, while requiring that their asymptotic fits smoothly join the 
the app and a vv total cross sections at ^fs =4.0 GeV — again both in magnitude and slope. In both cases, 
we have used all of the extensive data of the PDG group [TU). However, we then subject these data to a 
screening process, the "Sieve" algorithm^, in order to eliminate "outliers" that can skew a \ 2 fit- With 
the "Sieve" algorithm, a robust fit using a Lorentzian distribution is first made to all of the data to sieve 
out abnormally high A% 2 , the individual i th point's contribution to the total \ 2 ■ The \ 2 fits are then 
made to the sieved data. Both the ixp and nucleon-nucleon systems strongly favor a high energy In 2 s 
fit of the form: a ± = c + ciln (^) + c 2 ln 2 (^) + /3 P , ± S basically excluding a Ins 

fit of the form: a ± = Co + ciln (^A + /3-pi 1 ± <5 (t^)" • The upper sign is for ir + p (pp) and the 

lower sign is for ir~p (pp) scattering, where v is the laboratory pion (proton) energy, and m is the pion 
(proton) mass. 



High energy cross sections for the scattering of hadrons should be bounded by a ~ In 2 s, where s is the 
square of the cms energy. This fundamental result is derived from unitarity and analyticity by Froissart 2 , 
who states: "At forward or backward angles, the modulus of the amplitude behaves at most like sin 2 s, as 
s goes to infinity. We can use the optical theorem to derive that the total cross sections behave at most like 
In 2 s, as s goes to infinity". In this context, saturating the Froissart bound refers to an energy dependence 
of the total cross section rising no more rapidly than In 2 s. 

The question as to whether any of the present day high energy data for pp, pp and tt + p, tt~p cross 
sections saturate the Froissart bound has not been settled; one can not unambiguously discriminate between 
asymptotic fits of Ins and In 2 s using high energy data only [20]. We here point out that this ambiguity 
is resolved by requiring that the fits to the high energy data smoothly join the cross section and energy 
dependence obtained by averaging the resonances at low energy. Imposing this duality 5 condition, we show 
that only fits to the high energy data behaving as In 2 s that smoothly join (in both magnitude and first 
derivative) to the low energy data at the "transition energy" (defined as the energy region just after the 
resonance regions end) can adequately describe the highest energy points. This technique has recently been 
successfully used by Block and HalzenjH] to show that the Froissart bound is saturated for the jp system. 

We will use real analytic amplitudes to describe the data. Following Block and CahnjJJ, we write the 
crossing-even real analytic amplitude for high energy scattering asj^jEl 



/+ = ^{ A + PIMs/so) - W2] 2 + af-ieMi-ri/* _ »^/ + (0)| , 



(1) 



and the crossing-odd amplitude as 

/_ = -Da a_1 e* ff(1-a/2 . (2) 

where A, fj, c, sq and p are real constants. The variable s is the square of the center of mass system (cms) 
energy and v is the laboratory momentum. The additional real constant /+(0) is the subtraction constant 
at v — needed to be introduced in a singly-subtracted dispersion relational , [SI ■ Using the optical theorem, 
we obtain the total cross section 



a ± = A + 13 



In s/s - — 



c sm(-Kp/2)s^ 1 ± D cos(7ra/2)s a_1 (3) 



with p, the ratio of the real to the imaginary part of the forward scattering amplitude, given by 

p ± = — j/3 7rlns/s -ccos(7r^/2)s^ 1 + — / + (0)±Z?sin( 7 ra/2)s Q - 1 l, (4) 

Otot I V ) 

where the upper sign is for ir + p (pp) and the lower sign is for ir~p (pp) scattering, and the even amplitude 
applies to the spin-averaged jp scattering 

We now introduce the definitions A = cq + — j^, so = 2m 2 e~ Cl ^ 2c2 \ (3 — C2, c = ^7^7^ Pv and 
D = ( 2m , ^ 5. In the high energy limit, where s — > 2mv, Eq. © and Eq. along with their cross section 



derivatives -jr 2 -, — r, can be written as 



a ± 


= co-f 




1 


± 




'm) 





Cl ln(^) + c 2 ln*(^) + fo,(^y-±6(^-Y-\ (5) 
Vm/ \m/ \m/ Vm/ 

- Cl+C27 rln(-) -/Sp/cot (-f (- +— /+ )±*tan — (- , 6 

2 \mJ 2 \mJ v 2 \m/ J 

da ± ( 1 1 f21n(^/m))l ., , 91 . , 

J 1 '-02 < +/v ^-i)wr 2 (7) 



(v/m)) J (v/m)) 

±5{(a-l)(y/m)r- 3 }, (8) 

where the upper sign is for n + p (pp) and the lower sign is for Tr~p (pp) scattering. The exponents p and 
a are real. This transformation linearizes Eq. (J5J in the real coefficients cq,ci,C2,0p' an d S, convenient for 
a x 2 n t to the experimental total cross sections and p-values. Throughout we will use units of v and m in 
GeV and cross section in mb, where m is the projectile mass. 
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Let <t + be the total cross section for it + p (pp) scattering and a the total cross section for 7r p (pp) 
scattering. It is convenient to define, at the transition energy isq, 

a + {v /m) + a~(v /m) 
<x av = 

= c + ci ln(v /m) + c 2 \n 2 (u /m) + /V (t'o/m)^ 1 , (9) 



Act 



777. a 



<7 + (v /m) - a (i/ /m) 



2 

\Q-1 



0(7/0/™)" , (10) 
1 / do~ + do~~ 



A? 



2 ^d(i//m) d{v/m)J v=vq 

1 1 | f 21n(m/m) 

(vo/m) J 2 \ (yo/m)) 

1 / G?CT + d(T _ 



Cl 



| + / 3p,{( M -l)(Vm)' 1 - 2 }, (11) 



2 \d(v/m) d(v/m)J v=Vq 
= 6{(a-l)(v a /m) a - 2 }. (12) 

Using the definitions of cr av , Ac, 777 av and Am, we now write the four constraint equations 



(Vrn) 2 -f 

P7" = : 

fl- 1 



1 1 (2\xi{v Q /m) 

m av - Cl < 7 ; r > - C 2 



(1/0/m) J ' 1 {va/m) 



(13) 



co = cr av - ci \n(v Q /m) - c 2 In (vo/m) - (3 v >{v /my , (14) 

a = l + -r-(Mo/m), (15) 
Act 

6 = ACTS/™) 1 "", (16) 

that utilize the two slopes and the two intercepts at the transition energy v>q, where we join on to the 
asymptotic fit. We pick Vq as the (very low) energy just after which resonance behavior finishes. We use 
l-i = 0.5 throughout, which is appropriate for a Regge-descending trajectory. In the above, m = m p is the 
proton mass for the pp and pp systems, while m = m^ is the pion mass for the ir~p and n + p systems. 

Our strategy is to use the rich amount of low energy data to constrain our high energy fit. At the 
transition energy i^o, the cross sections a + (vQ/m) and er~(7/o/m), along with the slopes ^ d(J/m) ) ana - 

( d{v/m) ) ' are use d 1° constrain the asymptotic high energy fit so that it matches the low energy data 
at the transition energy v$. We pick v$ much below the energy at which we start our high energy fit, but 
at an energy safely above the resonance regions. Very local fits are made to the region about the energy 
7/0 in order to evaluate the two cross sections and their two derivatives at vq that are needed in the above 
constraint equations. We next impose the 4 constraint equations, Equations (|13fl . (|14fl . (|15fl and (|16fl . which 
we use in our \ 2 fit to Equations [5] and [BJ For safety, we start the data fitting at an energy z/ m i n appreciably 
higher than the transition energy. The transition energies, with appropriate cross sections and slopes, are 
summarized in Tabic ^ along with the minimum energies used in the asymptotic fits. 

We stress that the odd amplitude parameters a and S and hence the odd amplitude itself is completely 
determined by the experimental values Am and Act at the transition energy vq. Thus, at all energies, the 
differences of the cross sections ct~ — ct + (from the optical theorem, the differences in the imaginary portion 
of the scattering amplitude) and the differences of the real portion of the scattering amplitude are completely 
fixed before we make our fit. Further, for a In 2 s (Ins) fit, the even amplitude parameters Co and f3' v are 
determined by c\ and C2 (ci only) along with the experimental values of CT av and m av at the transition energy 
vq. In particular, for a In 2 s (Ins) fit, we only fit the 3 (2) parameters c±, C2, and /(0) (ci and /+(0)). Since 
the subtraction constant /+(0) only enters into the p- value determinations, only the 2 parameters c\ and Ci 
of the original 7 are required for a In 2 s fit to the cross sections ct^ , which gives us exceedingly little freedom 
in this fit — it is indeed very tightly constrained, with not much latitude for adjustment. The cross sections 
ct for the Ins fit are even more tightly constrained, with only one adjustable parameter, c\. 
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We now outline the adaptive Sieve algorithm^ that minimizes the effect that "outliers" — points with 
abnormally high contributions to x 2 — have on a fit when they contaminate a data sample that is otherwise 
Gaussianly distributed. Our fitting procedure consists of several steps: 

1. Make a robust fit of all of the data (presumed outliers and all) by minimizing Aq, the Lorentzian 
squared with respect to a, where 

N 

A 2 a ( a ;x) = £ln{l + 0.179A X ?(a: i ;a)} , (17) 

with a = {aj., . . . , ojf } being the M-dimensional parameter space of the fit. x — {x\, . . . ,Xn} 
represents the abscissa of the N experimental measurements y = {yi, . . . , yjv} that are being fit and 

A-Xi( x i'i a ) = ^ y '~ v ^' ,a ^ is the individual x 2 contribution of the i th point, where y(xf,a) is the 
theoretical value at Xi and Ui is the experimental error. It is shown in ref. yQ that for Gaussianly 
distributed data, minimizing Aq gives, on average, the same total x^in = Ei=i ^X 2 i x u a ) from Eq. ((T7|) 
as that found in a conventional x 2 fit, as well as rms widths (errors) for the parameters that are almost 
the same as those found in a \ 2 fit- 

A quantitative measure of whether point i is an outlier, i.e., whether it is "far away" from the true 

signal, is the magnitude of its Axf (2^; a) = ^ t/ '~ y ^" ct ^ . The reason for minimizing the Lorentzian 

squared is that this procedure gives the outliers much less weight w in the fit (w oc 1/ \J Axf (xf, a), for 
large Axf (xn a)) than does a \ 2 fit {w oc yAxf ( x i\ a ))i thus making the fitted parameters insensitive 
to outliers and hence robust. For details, see ref. [Q. 

If Xmin i s satisfactory, make a conventional x 2 fit to get the errors and you are finished. If Xmin i s n °t 
satisfactory, proceed to step |21 

2. Using the above robust Aq fit as the initial estimator for the theoretical curve, evaluate Axf(£,;a), 
for each of the N experimental points. 

3. A largest cut, Axf {xi\ a) max , must now be selected. We start the process with Axf «) max = 9- If 
any of the points have Ax 2 (xi]a) > Axf (xi\ a) max , reject them — they fell through the "Sieve". The 
choice of Axf (xf, a) max is an attempt to pick the largest "Sieve" size (largest Ax 2 (xi; a) max ) that 
rejects all of the outliers, while minimizing the number of signal points rejected. 

4. Next, make a conventional x 2 fit to the sifted set — these data points are the ones that have been 
retained in the "Sieve". This fit is used to estimate Xmin- Since the data set has been truncated by 
eliminating the points with Ax 2 (xi] a) > Axf {xi\ a) max , we must slightly renormalize the Xmin found 
to take this into account, by the factor 1Z. For Axf max = 9, 6, and 4, the factor 1Z is given by 1.027, 
1.140 and 1.291, whereas the fraction of the points that should survive this x 2 cu t — for a Gaussian 
distribution — is 0.9973, 0.9857 and 0.9545, respectively. A plot of TZ" 1 as a function of Axf max is given 
in Figure ^ which is taken from ref. £Q. 

If the renormalized Xmin' *- e -> ^ x Xmin i s acceptable — in the conventional sense, using the ordinary 
X 2 distribution probability function — we consider the fit of the data to the model to be satisfactory 
and proceed to the next step. If the renormalized Xmin ^ s n °t acceptable and Axf (xf, a) m ax is not too 
small, we pick a smaller Axf {xi\ a) m ax and go back to step [3] The smallest value of Axf (xf, a)max 
that we used is Axf (x^, a) max = 4. 

5. From the x 2 fit that was made to the "sifted" data in the preceding step, evaluate the parameters 
a. Next, evaluate the M x M covariance (squared error) matrix of the parameter space which was 
found in the x 2 fit- We find the new squared error matrix for the A 2 fit by multiplying the covariance 
matrix by the square of the factor r x 2. From Figure we find that r x 2 ~ 1.02,1.05 and 1.11 for 
Axf (xf, a)max = 9, 6 and 4, respectively . The values of r x 2 > 1 reflect the fact that a x 2 fit to the 
truncated Gaussian distribution that we obtain — after first making a robust fit — has a rms (root mean 
square) width which is somewhat greater than the rms width of the x 2 fit to the same untruncated 
distribution 1 . 
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The application of a x 2 fit to the sifted set gives stable estimates of the model parameters ex, as well as a 
goodncss-of-fit of the data to the model when Xmin i s renormalized for the effect of truncation due to the cut 
A\ 2 (xi; a) max . One can now use conventional probabilities for x 2 fits, i-e., the probability that \ 2 is greater 
than 1Z x Xmini f° r the number of degrees of freedom v. Model parameter errors are found by multiplying 
the covariance (squared error) matrix of the conventional x 2 fit by the appropriate factor {r x ^) 2 for the cut 
Axf (xi;a) max . 

Table [3 summarizes the results of our simultaneous fits to all of the available data from the Particle Data 
Group 10 for ov+p, cr^- p , p^+p and p v - p , using the 4 constraint equations with a transition energy y/s = 2.6 
GeV and a minimum fitting energy of 6 GeV, after applying the "Sieve" algorithm^. Three Ax 2 max cuts, 
4, 6 and 9, were made for Xv^^/m^) fits. There was considerable improvement in the renormalized x 2 /d.f. 
going from Ax 2 max = 9 to Ax 2 max = 6. However, there was no improvement of the renormalized x 2 /d.f. 
going from Ay? D „ = 6 to AY- mlv = 4 — indeed, it increased from 1.294 to 1.364. Since the errors also 
become substantially larger for the Ax 2 max = 4 cut, we chose to use the values of the In (u/m^) fit with a 
Ax 2 max = 6 cut. This cut was therefore also used for the ln(i//m w ) fit. The probability of the fit for the 
data set using the Ax 2 max = 6 cut was ~ 0.02, a somewhat low probability, albeit one that is often deemed 
acceptable in a fit with this many degrees of freedom (d.f.=127). In contrast, the probability of the hiiy/m^) 
fit using the Ax 2 max = 6 data set is << 10 -16 and is clearly ruled out, as is graphically demonstrated in 
Fig.H 

It should be noted that when using a ln 2 (i//m„-) fit before imposing the "Sieve" algorithm, a value of 
X 2 /d.f.=3.472 for 152 degrees of freedom was found, compared to x 2 /d.f.=1.294 for 127 degrees of freedom 
when using the A%f max = 6 cut. In essence, the "Sieve" algorithm eliminated 25 points with energies y/s > 6 
GeV (2 cr w + p , 19 o v - v , 4 p n + p ), while changing the total renormalized x 2 from 527.8 to 164.3. These 25 
points that were screened out had a x 2 contribution of 363.5, an average value of 14.5. If the distribution 
had been Gaussian with no outliers, one would have expected about 2 points having Axf > 6, giving a total 
X 2 contribution slightly larger than 12, compared to the observed value of 363.5. Thus, we see the effect of 
the "Sieve" algorithm in cleaning up the data sample by eliminating the outliers. 

Next, we analyze the pp and pp systems. Table |3 summarizes the results of our simultaneous fits to the 
available accelerator data from the Particle Data Groups]] for <j pp , o pp , p pp and p pp , using the 4 constraint 
equations with a transition energy y/s = 4 GeV and a minimum fitting energy of 6 GeV, again using the 
"Sieve" algorithm. Two Ax 2 max cuts, 6 and 9, were made for ln 2 (^/m p ) fits. The probability of the fit for 
the cut Ax 2 max = 6 was ~ 0.2, a very satisfactory probability for this many degrees of freedom, and we 
chose this data set rather than the data set corresponding to the Axf max = 9 cut. As seen in Table the 
fitted parameters are very insensitive to this choice. The same data set (Ax 2 max = 6 cut) was also used for 
the ln(^/m p ) fit. The probability of the ln(V/m p ) fit is << 10 -16 and is clearly ruled out. This is illustrated 
graphically in Fig.[|)]and Fig.0 

We note that when using a m 2 (V/m p ) fit before imposing the "Sieve" algorithm, a value of x 2 /d.f.=5.657 
for 209 degrees of freedom was found, compared to x 2 /d.f.=1.095 for 184 degrees of freedom when using the 
Ax 2 max = 6 cut. The "Sieve" algorithm eliminated 25 points with energies y/s > 6 GeV (5 a pp , 5 crp p , 15 
p pp ), while changing the total renormalized x 2 from 1182.3 to 201.4. These 25 points that were screened out 
had a x 2 contribution of 980.9, an average value of 39.2. For a Gaussian distribution, about 3 points with 
Ax? > 6 are expected, with a total x 2 contribution of slightly more than 18 and not 980.9. Again, we see 
the effect of the "Sieve" algorithm in ridding the data sample of outliers. 

Fig. [3shows the individual fitted cross sections (in mb) for n + p and 7r~p for In 2 [y / m ff ) and ln(i//m T ), 
for the cut Axf max = 6, from Table [21 plotted against the cms energy, y/s, in GeV. The data shown are 
the sieved data with y/s > 6 GeV. The In (v/m-,/) fits with Ax 2 max = 6 corresponding to the solid curve 
for n~p and the dash-dotted curve for ir + p, are in excellent agreement with the cross section data. On the 
other hand, the ln(i//m w ) fits — the long dashed curve for ir~p and the short dashed curve for 7r + p — although 
they fit in the low energy region almost identically to the ln 2 (j//m 7r ) fits — are very bad fits which clearly 
underestimate all of the high energy cross sections, leading to huge Xmirp an d hence are ruled out. 

Figure |3 shows the individual fitted p-values for ir + p and ir~p for ln 2 (^/m,r) and ln(i>/m v ), for the cut 
Ax 2 max = 6, from Tabled vs. y/s, the cms energy in GeV. The data shown are the sieved data with y/s > 6 
GeV. The ln 2 (i//m 7r ) fits with Axf max = 6, corresponding to the solid curve for ir~p and the dash-dotted 
curve for Tt + p, reproduce the data reasonably well. On the other hand, the ln(j//m w ) fits are rather poor. 
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Again the p data offer firm support for the Froissart bound fits, while ruling out ln(i//m n ) fits for the np 
system. 

Figure 0] shows an expanded scale of energies, in which all available 7175 cross sections are shown, from 
threshold to the highest available energies. The dashed curve is the even amplitude pion cross section 

y "" +p "^ T computed from our Axf max = 6 cut, whereas the solid curve is the result of a similar analysis 

for the spin-averaged (even) cross sectionjoj for jp scattering, rescaled by multiplying it by 2f0, a familiar 
number from the vector dominance model. It is most striking that these two independent curves are virtually 
indistinguishable in the entire energy interval in which experimental data are available, i.e., 2 < yfs < 300 
GeV — a result most strongly supporting the vector dominance model. 

All known np cross section data are plotted in Figure^ which compares our analysis using the 4 constraint 
equations (Axf max = 6 from Table EJl with the analysis of Igi and Ishida[S] which used finite energy sum 
rules (FESR) for their low energy data. They only fitted the even cross section, so we have plotted in Fig. 
Ufa) the even portion of our hi 2 (i>/m 7r ) fit as the solid curve. It is seen to go smoothly through the average 

cross section, { ^ 7 ' +p \ C ' 1 ' p ^j, for pion- proton scattering. The dashed-dot curve, using the FESR, is from Igi 

and Ishida[S]. It does not go very smoothly through the average of the points, but rather goes much closer 
to ov+p in the energy region from 10 to 30 GeV. Perhaps this is the result of their trying to fit only the even 
cross section, whereas we separately fit o~ w + p and av-p. We have plotted in Fig. |Sfb) the even portion of our 
ln(i//m T ) fit as the dashed curve, with the FESR result being the dashed-dot-dot curve. Clearly both curves 
rule out a ln(i//m T ) behavior. Both analyses strongly support a hi 2 (i> / m^) behavior and thus a saturation 
of the Froissart bound for the irp system. 

Figure shows the individual fitted cross sections (in mb) for pp and pp for ln 2 (z^/m p ) and ln(^/m p ) for 
the cut Axf max = 6 in Table |3J plotted against the cms energy, s/s, in GeV. The data shown are the sieved 
data with y/s > 6 GeV. The In 2 (v/m p ) fits to the data sample with Ax 2 max = 6, corresponding to the solid 
curve for pp and the dash-dotted curve for pp, are excellent, yielding a total renormalized x 2 — 201.5, for 184 
degrees of freedom, corresponding to a fit probability of ~ 0.2. On the other hand, the ln(v/m p ) fits to the 
same data sample — the long dashed curve for pp and the short dashed curve for pp — are very bad fits, yielding 
a total x 2 = 2613.7 for 185 degrees of freedom, corresponding to a fit probability of << 10~ 16 . In essence, 
the ln(^/m p ) fit clearly undershoots all of the high energy cross sections. The ability of nucleon-nucleon 
scattering to distinguish cleanly between an energy dependence of ln 2 (^/m p ) and an energy dependence of 
ln(^/m p ) is even more dramatic than the pion result. 

Figure {7\ shows the individual fitted p-values for pp and pp In 2 (y / m p ) and hx{v/m p ) from Table|3J using 
Axf max = 6 — plotted against the cms energy, s/s, in GeV. The data shown are the sieved data with y^i > 6 
GeV. The In 2 (y/m p ) fits, corresponding to the solid curve for pp and the dash-dotted curve for pp, fit the 
data reasonably well. On the other hand, the ln(u/m p ) fits, the long dashed curve for pp and the short 
dashed curve for pp, are very poor fits, missing completely the precise p pp at 546 GeV, as well as p pp at 1800 
GeV. These results again strongly support the \n 2 (u/m p ) fits that saturate the Froissart bound and once 
again rule out ln(i//m p ) fits for the pp and pp system. 

A few remarks on our ln 2 (^/m p ) asymptotic energy analysis for pp and pp are in order. It should be 
stressed that we used both the CDF and E710/E811 high energy experimental cross sections at \fs = 1800 
GeV in the ln 2 (i//m p ) analysis, summarized in Tabled Ax| max = 6 and shown in Figures|Uland[7| Inspection 
of Fig. shows that at y/s = 1800 GeV, our fit effectively passes below the cross section point of ~ 80 mb 
(CDF collaboration). In particular, to test the sensitivity of our fit to the differences between the highest 
energy accelerator pp cross sections from the Tevatron, we next omitted completely the CDF (~ 80 mb) 
point and refitted the data without it. This fit, also using Ax 2 max = 6, had a renormalized x 2 /d.f.=1.055, 
compared to 1.095 with the CDF point included. Since you only expect, on average, a Ax 2 of ~ 1 for the 
removal of one point, the removal of the CDF point slightly improved the goodness-of-fit. Moreover, the 
new parameters of the fit were only very minimally changed. As an example, the predicted value from the 
new fit for the cross section at = 1800 GeV — without the CDF point — was tr pp = 75.1 ± 0.6 mb, where 
the error is the statistical error due to the errors in the fitted parameters. Conversely, the predicted value 
from Table El— which used both the CDF and the E710/E811 point — was <7 pp = 75.2 ± 0.6 mb, virtually 
identical. Further, at y/s = 14 TcV (LHC energy), the fit without the CDF point had a pp = 107.2 ± 1.2, 
whereas including the CDF point (Table gave cr pp = 107. 3± 1.2. Thus, within errors, there was practically 
no effect of either including or excluding the CDF point. The fit was determined almost exclusively by the 
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E710/E811 cross section — presumably because the asymptotic fit was locked into the low energy transition 
energy vq , thus sampling the rich amount of lower energy data. 

Our result concerning the (un)importance of the CDF point relative to E710/E811 result is to be con- 
trasted with the statement from the COMPETE Collaboration^) which emphasized that there is: "the 
systematic uncertainty coming from the discrepancy between different FNAL measurements of <7tot" , which 
contribute large differences to their fit predictions at high energy, depending on which data set they use. 
In marked contrast to our results, they conclude that their fitting techniques favor the CDF point. Our 
results indicate that both the cross section and p-value of the E710/E811 groups are slightly favored. More 
importantly, we find virtually no sensitivity to high energy predictions when we do not use the CDF point 
and only use the E710/E811 measurements. Our method of fitting the data — by anchoring the asymptotic fit 
at the low transition energy v — shows that our high energy predictions are quasi-independent of the FNAL 
"discrepancy", leading us to believe that our high energy cross section predictions at both the LHC and 
at cosmic ray energies are both robust and accurate. In Table 0] we give predictions — from our ln 2 (i//m p ) 
fit — for some values of Op p and p pp at high energies. The errors quoted are due to the statistical errors of 
the fitted parameters ci, C2 and /+(0) given in the A X f max = 6, \n 2 {v/m p ) fit of Tabled! 

In Fig. |H1 we show an extended energy scale, from threshold up to cosmic ray energies (1.876 < y^s < 10 5 
GeV), plotting all available pp and pp cross sections, including cosmic ray pp cross sections inferred from 
cosmic ray p-air experiments by Block, Halzen and Stanov 11 . The solid curve is our result from Table [3] 
of the even cross section from \x?(yjm p \ A% 2 = 6. The dashed-dot-dot curve is from an independent 
QCD-inspired eikonal analysis |llj of the nucleon-nucleon system. The agreement is quite remarkable — the 
two independent curves are virtually indistinguishable over almost 5 decades of cms energy, from ~ 3 GeV 
to 100 TeV. Figure [S] clearly indicates that the pp and pp cross section data greater than ~ 3 GeV can be 

explained by a fit of the form — cq + c\ In ( — \ + C2 In 2 ( — \ + B-p* { — \ ± 6 ( — \ over an 

\ TTLp J \ TTLp J I ITLp J \ TTlp J 

enormous energy range, i.e., by a In 2 s saturation of the Froissart bound. 

In Table we make predictions of total cross sections and p- values for pp and pp scattering — in the low 
energy regions covered by RHIC, together with the energies of the Tevatron and LHC as well as the high 
energy regions appropriate to cosmic ray air shower experiments. 

We give strong support to vector meson dominance by showing that the even cross section from our fits 
for ir + p and ir~p data agrees exceedingly well with a rescaled (multiplied by a factor of 210) a lp analysis 
done earlier by Block and Halzen[5|, when both cross sections have a In 2 s asymptotic behavior. 

In conclusion, we have demonstrated that the duality requirement that high energy cross sections 
smoothly interpolate into the resonance region strongly favors a In 2 s behavior of the asymptotic cross sec- 
tions for both the irp and nucleon-nucleon systems, in agreement with our earlier result for jp scattering 6 . 
We conclude that the three hadronic systems, jp, irp and nucleon-nucleon, all have an asymptotic In 2 s 
behavior, thus saturating the Froissart bound. 

At 14 TeV, we predict Cpp — 107.3zbl.l mb and ppp — 0.132±0.001 for the Large Hadron Collider— robust 
predictions that rely critically on the saturation of the Froissart bound. 
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Transition Energy Parameters 


7r + p and 7r p Scattering 


pp and pp Scattering 


vo, lab transition energy (GeV) 


3.12 


7.59 


— > \/s , cms transition energy (GeV) 


2.6 


4 


(J + (v ) (mb) 


28.91 


40.18 




32.04 


56.99 


UCM))^, (mb) 


-0.2305 
-1.446 


-0.2262 
-0.2740 


Minimum fitting energy 




^ m i n , lab minimum energy (GeV) 


18.71 


18.25 


— > Vs m in' cms minimum energy (GeV) 


6.0 


6.0 



m is the pion (proton) mass and v is the laboratory pion (proton) energy 



Table 1: The transition energy parameters used for fitting 7r + p, tv p , pp and pp scattering. 



8 





(J <~ ^(^/TOtt) 


a ~ \n(v/m^) 


Parameters 


Ay 2 

Az max 


Ay 2 

a. i max 




6 


9 


6 


Even Amplitude 


c (mb) 


20.11 


20.32 


12.75 


ci (mb) 


-0.921 ±0.110 


-0.981 ±0.100 


1.286 ±0.0056 


c 2 (mb) 


0.1767 ±0.0085 


0.1815 ±0.0077 




f3 r , (mb) 


54.40 


54.10 


64.87 


M 


0.5 


0.5 


0.5 


/(0) (mb GcV) 


-2.33 ±0.36 


-2.31 ±0.35 


0.34 ±0.36 


Odd Amplitude 


5 (mb) 


-4.51 


-4.51 


-4.51 


a 


0.660 


0.660 


0.660 


x 2 - 

A mm 


148.1 


204.4 


941.8 


^ X Xmin 


164.3 


210.0 


1044.9 


V (d.f). 


127 


135 


128 


^ X Xmk> 


1.294 


1.555 


8.163 



Table 2: The fitted results for a 3-parameter \ 2 fit with a ~ ln 2 (^/m w ) and a 2-parameter fit with a ~ ln(i//m^-) 
to the total cross sections and p- values for 7r + p and 7r~p scattering. The renormalized xV^mm, taking into account 
the effects of the Ax? max cut, is given in the row labeled 1Z x Xmin/^- The errors in the fitted parameters have been 
multiplied by the appropriate r X 2. The pion mass is and the laboratory pion energy is v. 
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a ~ In (f/m p ) 


a ~ \sx{y/m p ) 


Parameters 


Ax; 


max 


Ay 2 

A. i max 




6 


9 


6 


Even Amplitude 


c (mb) 


37.32 


37.25 


28.26 


ci (mb) 


-1.440 ±0.070 


-1.416 ±0.066 


2.651 ±0.0070 


c 2 (mb) 


0.2817 ±0.0064 


0.2792 ±0.0059 




/V (mb) 


37.10 


37.17 


47.98 


M 


0.5 


0.5 


0.5 


/(0) (mb GeV) 


-0.075 ±0.59 


-0.069 ±0.57 


4.28 ±0.59 


Odd Amplitude 


8 (mb) 


-28.56 


-28.56 


-28.56 


a 


0.415 


0.415 


0.415 


x 2 - 

A mm 


181.6 


216.6 


2355.7 


^ X Xmin 


201.5 


222.5 


2613.7 


V (di). 


184 


189 


185 




1.095 


1.178 


14.13 



Table 3: The fitted results for a 3-parameter % 2 fit with a ~ ln 2 (^/m p ) and a 2-parameter fit with cr ~ m(z^/m p ) 
to the total cross sections and p- values for pp and pp scattering. The renormalized X^/i/nun, taking into account the 
effects of the Axi max cut, is given in the row labeled 1Z x Xmin/^- The errors in the fitted parameters have been 
multiplied by the appropriate r X 2. The proton mass is m p and the laboratory nucleon energy is v. 



,/s, in GeV 


(jpp, in mb 


Ppp 


a pv , in mb 


Ppp 


6 


48.97 ±0.01 


-0.087 ±0.008 


38.91 ±0.01 


-.307 ±0.001 


60 


43.86 ±0.04 


0.089 ±0.001 


43.20 ±0.04 


0.079 ±0.001 


100 


46.59 ±0.08 


0.108 ±0.001 


46.23 ±0.08 


0.103 ±0.001 


300 


55.03 ±0.21 


0.131 ±0.001 


54.93 ±0.21 


0.130 ±0.002 


400 


57.76 ±0.25 


0.134 ±0.002 


57.68 ±0.25 


0.133 ±0.002 


540 


60.81 ±0.29 


0.137 ±0.002 


60.76 ±0.29 


0.136 ±0.002 


1,800 


75.19 ±0.55 


0.139 ±0.001 


75.18 ±0.55 


0.139 ±0.001 


14,000 


107.3 ± 1.2 


0.132 ±0.001 


107.3 ± 1.2 


0.132 ±0.001 


16,000 


109.8 ± 1.3 


0.131 ±0.001 


109.8 ± 1.3 


0.131 ±0.001 


50,000 


132.1 ±1.7 


0.124 ±0.001 


132.1 ± 1.7 


0.124 ± 0.001 


100,000 


147.1 ±2.0 


0.120 ±0.001 


147.1 ±2.0 


0.120 ± 0.001 



Table 4: Predictions of high energy pp and pp total cross sections and p-values, from Table [3 a ~ ln 2 (z//m 7r ), 

max ^" 
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Figure 1: a) A plot of 1Z~ 1 , the reciprocal of the factor that multiplies X^ n [ n / U found in the x 2 fit to the sifted data set vs. 
Ax 2 cut, i.e., Ax 2 max - b) A plot of r x 2, the factor whose square multiplies the covariant matrix found in the x 2 fit to the 
sifted data set vs. Ax 2 cut, i.e., Ax 2 max - These figures are taken from ref. [Q. 
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10 100 

Figure 2: The fitted total cross sections cr^ +p and <r w - p in mb, vs. ^/s, in GeV, using the 4 constraints of Equations 1131 . 
1141 . 1151 and 1161 . The circles are the sieved data for n~p scattering and the squares are the sieved data for 7r + p scattering 
for > 6 GeV. The dash-dotted curve (n + p) and the solid curve (tt~p) are x 2 fits ("Table 121 a ~ ln 2 (^/m w ), A% 2 max = 6) of 

the high energy data of the form : cr 7v ± p = co + ciln (^) + c 2 ln 2 (^) + j3 v , (^-) M 1 ± <5 (^j)™ • The upper sign is for ir+p 
and the lower sign is for TT~p scattering. The short dashed curve (ir+p) and the long dashed curve (ir~p) are x 2 fits (Table 121 
a ~ ln(j//m w ), A X 2 max = 6 ) of the high energy data of the form : a 7T±p = c + Cl ln (^) + (^) M ~ 1 ± <5 ■ The 

upper sign is for 7r+p and the lower sign is for ir~p scattering. The laboratory energy of the pion is v and m is the pion mass. 
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Figure 3: The fitted p- values, p w + p and p^-„, vs. y/s, in GeV, using the 4 constraints of Equations 1131 . 1141 . 1151 and 
fl6t . The circles are the sieved data for it p scattering and the squares are the sieved data for 7r + p scattering for y/s > 6 
GeV. The dash-dotted curve (n + p) and the solid curve (tt~p) are \ 2 fits (Table 121 a ~ ln 2 (v/m w ), Ax^ m = 6) of the 

high energy data of the form : p± = ^ | f C1 + C2 tt In (^) - 0p, cot(7rp/2) (^) M_1 + ^/+(0) ± <5 tan(7ra/2) (£) a_1 }. 

The upper sign is for and the lower sign is for ir~ p scattering. The short dashed curve (tt+p) and the long 

dashed curve (ir~ p) are \ 2 fits (Table |21 a ~ ln(^/m^), Ax^ max = 6 ) of the high energy data of the form : p^ = 

j^lJ Cl ~ $P' cot(-7rp/2) + ~"/+(0) i <5tan(7ra/2) (~) j 1 - The upper sign is for w + p and the lower sign is 

for n~p scattering. The laboratory energy of the pion is v and m is the pion mass. 
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Figure 4: The circles are the cross section data for 7r p scattering and the squares are the cross section data for 7r+p 
scattering, in mb, vs. y/s, in GeV, for all of the known data. The dashed curve is the x 2 fit (Tabic a ~ ln 2 (v/m T ), 
Ax 2 max = 6) to the high energy cross section data of the even amplitude cross section, of the form : <T wp cvon = co + ciln (^j + 

C2ln 2 (^1 + fi-pi (7^)' i with co and ft-pi constrained by Eq. 1131 and Eq. 1141 , The laboratory energy of the pion is v 
and m is the pion mass. The dashed curve is 210 X cr JP , from a fit of jp cross sections by Block and Halzen J)j of the form: 
<t 7p = co + ciln(u/m p ) + C2m 2 (^/m p ) + fl-pi / W u/m p , where m p is the proton mass. The 7p cross sections were fit for cms 
energies y/s > 2.01 GeV, whereas the irp data (cross sections and p-values) were fit for cms energies y/s > 6 GeV. The two 
fitted curves are virtually indistinguishable in the energy region 2 < < 300 GeV. 
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Figure 5: The circles are the cross section data for it p scattering and the squares are the cross section data for n + p scattering 
for all known data, vs. y/s, in GeV. The solid curve in Fig. (a) is the \ 2 fit f Table |21 a ~ \v?(y /m^), Ax 2 max = 6) to the high 

energy cross section data of the even amplitude, of the form : eleven = cq + ciln (^) + C2ln 2 f^) + /Jp' (m)' , with co and 
/3-p/ constrained by Eq. 1131 and Eq. 1141 . The dash-dotted curve is an even amplitude In 2 (i//m n ) fit made by Igi and Ishida 5 , 
using finite energy sum rules (FESR). The dashed curve in Fig.(b) is the x 2 fit ( Tabled a ~ ln(u/m„), Ax 2 max = 6) to the 

high energy cross section data of the even amplitude, of the form : eleven = co + ciln (^)) + Pv ( ™) M > with Co an d P-pi 
constrained by Eq. 1131 and Eq. 1141 . The dot-dot-dashed curve is an even amplitude Inly/m^) fit made by Igi and Ishida 5 , 
using finite energy sum rules. The laboratory energy of the pion is v and m is the pion mass. 
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Figure 6: The fitted total cross sections a pp and up p in mb, vs. ^fs, in GeV, using the 4 constraints of Equations 1131 , 1141 . 
1151 and 1161 . The circles are the sieved data for pp scattering and the squares are the sieved data for pp scattering for y/s > 6 
GeV. The dash-dotted curve (pp) and the solid curve (pp) are \ 2 fits fTablel3l a ~ ln 2 (^/m,r), Ax 2 max = 6) of the high energy 

data of the form : cr^ = co + ci In (~) + C2ln 2 (~) + fi-p' (m) ^ -^(m) ' The upper sign is for pp and the lower sign is 
for pp scattering. The short dashed curve (pp) and the long dashed curve (pp) are \ 2 n * s ( Table l3l a ~ ln(^/m„-), Ax 2 max = 6 

) of the high energy data of the form : cr^ = cq + ciln + /3-pi (^) M i£ (m) ' u PP er s ig n is f° r PP an d the lower 
sign is for pp scattering. The laboratory energy of the nucleon is v and m is the nucleon mass. 
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Fi gUX6 7; The fitted /^-values, p pp and ppp. vs. y/s^ in GeV, using the 4 constraints of Equations 1131 . 1141 , 1151 and 
(Tjj). The circles are the sieved data for pp scattering and the squares are the sieved data for pp scattering for y/s > 6 
GeV. The dash-dotted curve (pp) and the solid curve (pp) are \ 2 fits fTablc 131 a ~ ln 2 (i'/m^-), Ax 2 max = 6) of the 

high energy data of the form : p± = | fci + c 2 7r In (i) - (i v , cot(7rp/2) {j£f~^ + £/+(0) ± $tan(7ra/2) (£) a_1 j. 

The upper sign is for pp and the lower sign is for pp scattering. The short dashed curve (pp) and the long dashed 
curve (pp) are x 2 fits (Table 151 a ~ ln^/m-n-), A%? max = 6 ) of the high energy data of the form : = 

^5 Cl — ^f" c °t( 7r / i /2) (^) M + ~77"/+(0) i <5tan(7ra/2) (^-) j 1 - The upper sign is for pp and the lower sign is for pp 
scattering. The laboratory energy of the nucleon is v and m is the nucleon mass. 
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^ s, in GeV 



Figure 8: The circles are the cross section data for pp scattering and the squares are the cross section data for pp scattering, in 
mb, vs. \/s, in GeV, for all of the known accelerator data. The solid curve is the \ 2 fit ((Tabled <r ~ \v? (v / m^) , Ax 2 max = 6) 

of the high energy data of the crossing-even amplitude, of the form : cr nn even = cq + ciln + C2ln 2 {^-^ + ftv 1 (m)^ ' 

with Co and /3-pi constrained by Eq. 1131 and Eq. 1141 . The dot-dot-dashed curve is the crossing-even amplitude cross section 
ann, from a QCD-inspircd fit that fit not only the accelerator pp and pp cross sections and p-values, but also fit the AGASA 
and Fly's Eye cosmic ray pp cross sections shown in the figure — work done several years ago by the Block, Halzen and Stanov 
(BHS group) EI. The laboratory energy of the proton is v and m is the proton mass. It is most striking that the two fitted 
curves for <r nll even, using on the one hand, the In 2 (u/m) model of this work and on the other hand, the QCD-inspired model of 
the BHS group lTTl . are virtually indistinguishable over 5 decades of cms energy, i.e., in the energy region 3 < y/s < 10 GeV. 
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